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Abstract Two different isospin splittings of nucleon effective mass in nuclear medium as the form of m? > m, and 
m < m, have been implemented in an isospin and momentum dependent transport model. Their impacts on the 
isospin emission in heavy-ion collisions is investigated thoroughly. It is found that the yield ratios of energetic 
neutrons to protons squeezed out during the compression stage of two colliding nuclides are sensitive to the isospin 
splitting. The elliptic flows of free nucleons are also to be promising observables for extracting the nucleon effective 
mass splitting. Further experimental measurements are being expected, in particular at the CSR-CEE platform in 
Lanzhou. Several observables are proposed for constraining the density dependence of symmetry energy, such as the 


transverse flow difference of neutrons and protons, double ratios of n/p and  / m`, excitation functions of z / 7° and 


KK’. 
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1 Introduction 


It has been well known that the masses of neutron and 
proton are equal each other in vacuum (about 1 GeV), 
and the effective mass in nuclear matter or finite nuclei 


é š L2 
deviates from its vacuum value!” 


Moreover, a 
splitting of neutron and proton effective mass exists in 
neutron-rich nuclear matter, which increases with the 
isospin asymmetry and nucleon density. Predictions of 
the mass splitting based on nuclear many-body 
differ 
nucleon-nucleon interactions, the Brueckner-Hartree- 
Fock (BHF) and Dirac-Brueckner-Hartree-Fock 


(DBHF) calculations predict a neutron-proton mass 
[3-6] 


theories widely. Based on the realistic 


splitting of m, > m, The models based on 
relativistic mean-field (RMF) theory give a contrary 
result after the inclusion of the isovector mesons p and 
8l, The Skyrme-Hartree-Fock (SHF) predicts both 
splittings of m, > m, and mÝ < m, exist with 
different Skyrme parameters!" More realistically, the 
mass splitting of neutron and proton has been 
constrained from the energy dependence of the Lane 


(symmetry) potential by the nucleon-nucleus 


scattering experimental data, which moderately 
supports the mass splitting of m,>m, M] Further 
constraints on the effective mass splitting from 
heavy-ion collisions are still necessary and its 
influence on reaction dynamics would be very 
interesting. 

The knowledge of the density dependence of 
the nuclear symmetry energy is still poorly known, in 
particular at high baryon densities, but which has 
important application in nuclear physics itself and also 
in astrophysics. Heavy-ion reactions with neutron-rich 
beams provide a unique opportunity to explore the 
density-dependent symmetry energy in a broad domain 
of density. In this work, we present systematic 
investigations of the effective mass splitting of neutron 
and proton and its influence on reaction dynamics. An 
extraction of the high-density symmetry energy from 
the preequilibrium nucleon emission and meson 
production (mand K ) is performed with an isospin and 


momentum dependent transport model!'!. 
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2 Model description 


The time evolutions of the baryons and mesons in the 


process of two colliding partners under the 
self-consistently generated mean-field are governed by 
Hamilton's equations of motion. The Hamiltonian of 
baryons consists of the relativistic energy, the effective 
interaction potential and the momentum dependent 
part!7! A Skyrme-type momentum-dependent nucleon 
-nucleon force distinguishing isospin effect is 
parameterized and implemented in the Lanzhou 
Quantum Molecular Dynamics (LQMD) model, which 
leads to a splitting of nucleon effective mass in 
neutron-rich matter. A Skyrme-type form for the 
momentum-dependent potential in the Hamiltonian by 
distinguishing isospin effect is taken in the LQMD 


model, which is expressed as 


lye, : 
Unon = 75 LC [[Japap'arf.(r,p) p 


f 2 f 
x| In(e(p—p’y +1)] frp’) 
Here C_ = (+x), Ca =Crom(l1 X) (TT) and 
the isospin symbols 7 ( T') represent proton or neutron. 


C mom mom ( 


The sign of x determines different mass splitting of 
proton and neutron in nuclear medium, e.g. positive 
signs corresponding to the case of mř < m’, . The 
parameters Cmom and € were determined by fitting the 
real part of optical potential as a function of incident 
energy from the proton-nucleus elastic scattering data. 
In the calculation, we take the values of 1.76 MeV, 
500 c?/GeV? for Cnom and s, respectively, which result 
in the effective mass m'/m=0.75 in nuclear medium at 
saturation density for symmetric nuclear matter. The 
parameter x can be adjusted as the strength of the mass 
splitting, and the values of —0.65 and 0.65 are 
respective to the cases of m,>m, and m,<m,. One 
can get the contribution of momentum-dependent 
interaction to symmetry energy for a cold nuclear with 
Fermi distribution for a phase-space density 
f.(r.p) = p,(r)@(p, (7) —[p))/ (4zp3(c)/3) with a 
Fermi momentum D,(T) = h(3x7p.)'°. The 
symmetry energy per nucleon in the LQMD model is 
composed of three parts, namely the kinetic energy, 
the local part and the momentum dependence of the 
potential energy as 


Eym(p)=55-( 3a e) + E55 (9) + Exe(p). (2) 


After the inclusion of momentum dependent 


contribution, a density, isospin and momentum 
dependent single-nucleon potential is obtained as 


follows: 


y 5/3 
8 oc 
U,(p,6,p)=a£ p| Z| +29, Z| ESOS 
p 3 LA 


0 Po 0 

, Esm (P) cco ôe? 
mpd + Ew (p)p 
op op, 


+ ma _Jev'f-@.p)[In(ep -p'? +] 
+ > [ot EpD[nep-p +] 
3) 


38l p4 8pp/ pP’, and d5°/dp,= 
—46p,/ p°. The nucleon effective (Landau) mass in 


Here TÆT', 


nuclear matter of isospin asymmetry 6=(Pn-Pp)/ P 
with on and pp being the neutron and proton density, 
respectively, is calculated through the potential as 


mm (1+7 


t 


with the free mass m, at 


Fermi momentum p=pr. Therefore, the nucleon 
effective mass only depends on the momentum- 
dependent term of the nucleon optical potential. 

Figure 1 is a comparison of the momentum 
dependence of single-nucleon optical potential with 
the mass splittings of m> m, in the left window and 
m, < m’, in the right window for the hard and 
supersoft symmetry energies, respectively. One should 
note that a cross of neutron and proton single-nucleon 
at high 


momentum, but which does not take place in the mass 


potentials appears in the case of m; > m’, 
splitting of m, <m’, and a broader separation exists 
with increasing the nucleon momentum. The 
difference of single-particle potentials for both mass 
splitting affects the isospin emission in heavy-ion 
collisions, in particular, in the domain of high-baryon 


densities. 
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k/ fm? 


k/ fm" 


Fig.1 (Color online)Momentum dependence of single-nucleon 
optical potential for isospin symmetric matter (d=0) and 
neutron-rich matter (60.2) with the mass splittings of 
m,>m, (left panel) and m,<m, (right panel), respectively. 

The scattering in two-particle collisions is 
performed by using a Monte Carlo procedure, in 
which the probability to be a channel in a collision is 
calculated by its contribution of the channel cross 
section to the total cross section. The primary products 
in nucleon-nucleon (NN) collisions in the region of 
1A GeV energies are the resonances of A(1232), 
N'(1440), N'(1535) and the pions. We have included 
the reaction channels as follows: 


NN eNA, NNe NN’, NN&AA, 
Ae Nz, N' & Nz, NN & NNz(s-state), (4) 
N’ (1535) > Np. 


At the considered energies, there are mostly A 
resonances which disintegrate into a z and a nucleon 
in the evolutions. However, the N” yet gives 
considerable contribution to the energetic pion yields. 

The strangeness is created as the secondary 
products in inelastic hadron-hadron collisions", We 
included the channels as follows: 


BB— BYK, BB—BBKK, Bz—- YK, 
Bz— NKK, Yr BK, BK >Yz, (5) 
YN — KNN. 


Here the B strands for (N, A, N^) and Y(A, È), K(K° ; 
K^) and K(K°, K). The elastic scattering between 
strangeness and baryons is considered through the 
channels of K BK B, Y B>Y Band K B> KB. The 
charge-exchange reactions between the K N-> K N and 


Y N>Y N channels are included by using the same 
cross sections with the elastic scattering, such as K? 
pK‘ n, K‘ n->K’ p ete. Correction of effective mass 
of kaons in nuclear medium on the elementary cross 
section is considered through the threshold energy, 
which results in the reduction of kaon and the 
enhancement of anti-kaon yields in heavy-ion 
collisions. 

The evolution of mesons (here mainly pions 
and kaons) is also determined by a Hamiltonian. Here, 
the Coulomb interaction is considered for pion and 
kaon propagation. The kaon and anti-kaon energy in 
the nuclear medium is computed from the chiral 
Lagrangian. A repulsive kaon-nucleon (KN) potential 
and attractive antikaon-nucleon potential with the 
values of 25.5 MeV and —96.8 MeV are used in the 
model, which result in the effective mass with 
m,./m,=1.05 and m,/m,=0.8 at saturation baryon 
density. 


3 Results and discussion 


The preequilibrium nucleons in high-energy heavy-ion 
collisions are mostly produced during a compression 
stage of two colliding partners within a very short time. 
Therefore, the high-density information of nuclear 
phase diagram is expected to be extracted from the 
nucleons or light complex particles, which can be 
constrained from the longitudinal rapidity distributions 
and the azimuthal emissions. Shown in Fig.2 is a 
comparison of transverse emission ratios of 
neutron/proton within the rapidity selection of 
\y/Ypro|<0.25 in the 197 Au+' Au reaction at the 
incident energy of 400A MeV for the near central 
(b=1 fm) and semi-central (b=6 fm) collisions, but 
with different mass splittings. One can see that the 
influence of the symmetry energy appears at low 
transverse momentum and the situation does not be 
changed with the mass splitting. A hard symmetry 
energy enforces a strongly repulsive force on neutrons 
in the high-density domain, furthermore, squeezes out 
more neutrons in the preequilibrium stage of 
dynamical evolution. From the negative contribution 
of the momentum-dependent interaction of m,>m, 
to the potential part of the symmetry energy in 
Ref.[12], one notices that the momentum-dependent 


potential leads to an attractive force on neutrons, in 
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particular at high densities, which reduces the n/p 
yields at high transverse momentum. Opposite 
contributions take place and a flat distribution appears 
in the case of m,<m,. The spectrum of n/p has a 
slightly different structure with the results in Ref.[15] 
that is caused from the momentum dependent 
contribution of symmetry potential. 
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Fig.2 Transverse momentum distributions of neutron/proton 
within the rapidity bin |y / Ypro;|<0.25 in the '”Au+'” Au reaction 
at the incident energy of 400 MeV/nucleon for the near central 
(b=1 fm) and semi-central (b=6 fm) collisions with different 
mass splitting. 

Collective flow has been verified as a nice 
approach to reconstruct the reaction plane and to study 
azimuthal correlation of the fireball formed in 
heavy-ion collisions. More sensitive observable can be 
seen from the flow difference between neutron and 
proton as shown in Fig.3. The m,> m, case is always 
larger, in particular in the domain of high transverse 
momentum. Basically, the symmetry energy does not 
change the spectrum structure. We calculated the 
difference of neutron and proton directed flows as a 
function of rapidity distribution in the '’Aut+'?’Au 
collisions (b=6 fm) at the energy of 400A MeV with 
the hard and supersoft symmetry energies, but for the 
different mass splitting as shown in Fig.4. One can see 
that an anti-flow spectrum appears in both mass 


splitting in comparison with nucleon flows!!617 


Moreover, the supersoft symmetry energy gives a 
larger flow difference, in particular around the region 
of projectile (target) rapidity. The phenomena can not 
be influenced by changing different mass splitting. The 
flows of fast nucleon emissions would be a nice probe 
for extracting the high-density behavior of nuclear 
symmetry energy. 
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Fig.3 Comparison of the difference between neutron and 
proton elliptic flows with the rapidity bin |y/Ypro|<0.25 for 
different mass splitting in semi-central '’Au+’’’Au collisions. 
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Fig.4 Transverse flow difference between neutrons and 
. . 197 197 “os 

protons in semi-central `° Au+ ‘Au collisions at the energy of 

400 MeV/nucleon with the hard and supersoft symmetry 

energies for different mass splittings of m,>m, and 

m; SIH 3 


Kaon meson is produced at the early stage in 
heavy-ion collision and promptly emitted after 
production, which can get directly the information of 
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high-density phase diagram. A pronounced effect of 
the stiffness of symmetry energy on kaon production 
can be observed from the spectrum of isospin ratio. 
The isospin effects appear at deep subthreshold 
energies as shown in Fig.5. The in-medium potential 
slightly changes the K°/K’ value because of its 
influence on the kaon propagation and also on the 
charge-exchange reactions. At the considered energies, 
the channel of NA->NYK contributes the main part for 
the kaon yields due to the larger production cross 
sections and the higher invariant energy, and the 
NN-NYK as well as zN YK have about one third 
contributions. One notices that a hard symmetry 
energy always has the larger values of the isospin 
ratios than the supersoft case in the domain of 
subthreshold energies (Ey,(K)=1.58 GeV). 


asy-hard (in-vacuum) 
asy-supersoft (in-vacuum) 
asy-hard (in-medium) 
asy-supersoft (in-medium) 


x 
x 
0.6 0.9 1.2 1.5 1.8 
E,,, / GeV-nucleon™* 
Fig.5 (Color online) Comparison of excitation functions of 


the K°/K" yields for central '’Au+'’Au collisions for the cases 
of hard and supersoft symmetry energies. 


4 Conclusion 


In summary, within the transport model (LQMD) we 
have investigated the impact of the momentum- 
dependent interaction on fast nucleon emissions in 
heavy-ion collisions. Two different mass splittings, i.e., 
mi > m, and mọ < m, , are chosen in the calculations. 
The momentum-dependent potential plays a 
significant role on the fast nucleon emissions. 
neutron/proton 


Specifically, the ratio at high 


transverse momenta is influenced, in which a flat 


spectrum appears in the case of m < m, . The elliptic 


flow difference between neutron and proton 
sensitively depends on the mass splitting, in particular 
at high transverse momenta. The neutron/proton ratio 
and the transverse flow difference of neutron and 
proton via transverse momentum distribution in the 
domain of mid-rapidity are sensitive to the stiffness of 
nuclear symmetry energy, which can be promising 
observables for extracting the high-density behavior of 
symmetry energy. The K°/K* ratio of neutron-rich 
heavy system in the domain of subthreshold energies 
is also sensitive to the stiffness of nuclear symmetry 


energy. 
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